Interplay between a hydrodynamic instability and a phase transition: the Faraday 

instability in dispersions of rodlike colloids 
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Strong effects of the Faraday instability on suspensions of rodlike colloidal particles are reported 
through measurements of the critical acceleration and of the surface wave amplitude. We show that 
the transition to parametrically excited surface waves displays discontinuous and hysteretic features. 
This subcritical behaviour is attributed to the shear-thinning properties of our colloidal suspensions 
thanks to a phenomenological model based on rheological data under large amplitude oscillatory 
shear. Birefringence measurements provide direct evidence that Faraday waves induce local nematic 
ordering of the rodlike colloids. While local alignment simply follows the surface oscillations for 
dilute, isotropic suspensions, permanent nematic patches are generated by surface waves in samples 
close to the isotropic-to-nematic transition and above the transition large domains align in the flow 
direction. This strong coupling between the fluid microstructure and a hydrodynamic instability is 
confirmed by numerical computations based on the microstructural response of rodlike viruses in 
shear flow. 

PACS numbers: 83.60.La, 83.50.Ax, 83.50.Rp 



INTRODUCTION 



Recently, the Faraday instability, i.e. the parametric 
instability of a fluid layer submitted to vertical vibra- 
tions [l], has appeared as a useful experiment to probe 
the interplay between the microstructure of various com- 
plex fluids and a classical hydrodynamic instability [2|- 
UM- Most previous experimental works have focused on 
the effects of viscoelasticity on the instability threshold 
0, HI, [l(|, on the subharmonic vs. harmonic response of 
the surface waves [H, EH , and on the influence of the mi- 
crostructure on the surface wave pattern [Sllllj. In shear- 
thinning semidilute solutions of polymers and of surfac- 
tant wormlike micelles, these effects are generally rather 
small perturbations to the classical instability of Newto- 
nian fluids. However, much stronger effects were found in 
a few other complex materials, namely clay suspensions 
[f|, shear-thickening cornstarch suspensions Q, and di- 
lute shear-thickening wormlike micelles [l2| , where local- 
ized hysteretic finger-like structures, persistent holes, and 
high-amplitude strip waves were respectively reported. 
These original observations were qualitatively interpreted 
0, @ and quantitatively modelled [l2} as a strong feed- 
back of the surface waves on the fluid microstructure: the 
oscillatory shear generated by small-amplitude perturba- 
tions triggers large variations of the shear rate-dependent 
viscosity that lead to discontinuous transitions to large- 
amplitude waves. 

Still, to the best of our knowledge, no direct exper- 
imental evidence for strong structural modifications in- 
duced by Faraday waves has been reported. The aim of 
the present study is to address this issue in suspensions of 
fd virus, a rodlike colloid that presents liquid crystalline 



phases |14Hl6l ]. These viruses are known to easily align 
under minute external perturbations and the local orga- 
nization of their microstructure can be simply probed by 
visualizing the birefringence field [I?], EH- Here, three 
suspensions of fd viruses at various concentrations, re- 
spectively deep into the isotropic domain, close to the 
isotropic-nematic transition, and in the nematic phase, 
are submitted to vertical vibrations. After a brief de- 
scription of the samples and of the experimental set up, 
we first report on surface wave amplitude measurements 
that present a strong hysteresis at the onset of the Fara- 
day instability. This subcritical behaviour is accounted 
for by a simple model based on the shear-thinning prop- 
erties of our rodlike colloids. We then move to bire- 
fringence experiments that reveal local nematic order- 
ing in vibrated suspensions that are initially isotropic 
at rest. These measurements are successfully compared 
to numerical simulations based on a Smoluchowski ap- 
proach 19|, 20[. Our results point to the existence of an 
out-of-equilibrium phase transition induced by a hydro- 
dynamic instability in suspensions of rodlike colloids. 



MATERIALS AND METHODS 

Our working fluids are suspensions of fd viruses in a 
buffer of 20 mM tris-HCl at pH=8.15 as described in 
Ref. [18j. The fd bacteriophage is a nearly perfectly 
monodisperse rodlike colloid of length 880 nm, diameter 
6.6 nm, and persistence length 2.2 /im. The number of el- 
ementary charges per unit length is around 1 e~.mii' 1 at 
pH=8.15. In order to screen the electrostatic interaction 
between viruses, NaCl is added to the buffer to obtain 
an ionic concentration of 20 mM. This system presents 
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FIG. 1: Schematic of the set up. A parallelepipedic cell (1) 
sealed by a plexiglass lid (2) and thermostated by water circu- 
lation (3) from a thermostated bath (4) is vertically vibrated 
by an electromagnetic shaker (5). The amplitude of the oscil- 
lation is recorded via an accelerometer (6), while the temper- 
ature of the cell is measured with a probe (7) and feedbacked 
to the water bath. Left: side view of the cell. 

various phases when the concentration c of fd virus is in- 
creased: isotropic, cholesteric, and smectic [ljl Eft Be- 
cause the free energy between the cholesteric phase and 
a nematic state is very low, we shall simply refer to the 
cholesteric phase as nematic in the following. We con- 
sider suspensions either in the isotropic or in the nematic 
phase by focusing on three different samples: an isotropic 
suspension of concentration c = 5.8 mg/mL well below 
the isotropic-nematic (I— N) phase transition, an isotropic 
suspension of concentration c = 11.3 mg/mL just be- 
low the I-N phase transition, and a nematic sample of 
concentration c = 13.6 mg/mL. When submitted to an 
external shear, fd viruses tend to align in the shear direc- 
tion. Deep into the isotropic phase, this alignment causes 
a dramatic decrease of the viscosity (shear-thinning be- 
haviour) |l7[ , while close to and above the I-N transition, 
more complex shear-induced phenomena such as vortic- 
ity banding, tumbling, and wagging have been reported 

m 

A dedicated set up (see fig. [J) was developed in order 
to measure both the birefringence intensity and the am- 
plitude of surface waves in a vertically vibrated layer of fd 
suspension. A parallelepipedic glass cell of length 72 mm, 
thickness 6 mm, and height 8 mm (respectively corre- 
sponding to the x, y, and z directions) is filled up to a 
height h = 6 mm with the suspension and vibrated by an 
electromagnetic shaker (Ling Dynamic Systems V406). 

The driving acceleration is sinusoidal with amplitude 
a = 10-100 m.s" 2 and frequency / = 50-200 Hz. The 
suspension is thermostated at T = 20 ± 0.5°C by a water 
circulation beneath the cell. The cell is sealed by a PVC 
cover to prevent evaporation and surface contamination. 
It is placed between two polarizers and lit from the back. 
A fast CCD camera (Mikrotron MC1310) captures the 
transmitted light at ~ 1000 fps. For each value of the 
acceleration a and frequency /, we wait for two minutes 
so that a steady state is reached, after which two movies 



are recorded over twenty driving periods (T = 20//): (i) 
a first movie with aligned polarizers from which the am- 
plitude t) is easily extracted as a function of the hori- 
zontal position x and time t and (ii) a second movie under 
crossed polarizers in order to measure the transmitted 
intensity field I(x, z, t) and to visualize the flow-induced 
birefringence.^ We checked that at the frequencies under 
study, the rather large aspect ratio of our experimental 
cell 72/6 > 10 always enforces a one-dimensional surface 
wave pattern along the x direction. In all cases, the sur- 
face response was observed to be subharmonic so that 
£(x, t) is well described by: 

t) = sin(fcE) ^2 £n cos ( y ( n+ ^) ujt + ' W 

where k is the wave number, u) = 1-nf, and (resp. ip n ) 
is the amplitude (resp. phase) of the n th harmonic. In 
our experiments £0 ^ £n>o , so that we will only consider 
n — in what follows. 

HYSTERESIS AND SUBCRITICALITY OF THE 
FARADAY WAVES 

Following classical protocols for detecting the onset of 
Faraday waves, we first determine the critical acceler- 
ation a c , i.e. the acceleration above which parametric 
waves first appear at the fluid surface at a given driv- 
ing frequency /, by (i) increasing the acceleration a by 
Sa/a = 0.1% every two minutes and noting the acceler- 
ation a" p when the surface first destabilizes, (ii) waiting 
for two more minutes for the surface to be fully desta- 
bilized and (iii) decreasing the acceleration by 8a every 
two minutes until the surface becomes flat again, which 
defines a second acceleration a^° wn . A striking difference 
with experiments in Newtonian fluids is the presence of 
a large hysteresis in the instability onset: in all three fd 
suspensions, a" p is about 4% greater than a^ own . For 
comparison, in a Newtonian fluid of similar viscosity, no 
hysteresis is detected, i.e. a" p = a^° wn to within our 
precision of 0.1%. 

The measurements of the surface wave amplitude £0 
shown in Fig. dja) for the more dilute sample reveal that 
this hysteresis is associated with a discontinuous transi- 
tion: above a" p , the wave amplitude jumps directly to 
several hundreds of micrometers and drops abruptly to 
zero when a is reduced by ~ 4%. As seen in Fig. [^b) 
and (c), hysteresis is also observed for the sample close 
to the I-N transition and for the nematic sample. More 
precisely, the hysteresis gets more pronounced as one gets 
closer to the I-N transition, and abruptly decreases after 
the transition. This evolution with c is reminiscent of 
that of the shear thinning properties of the sample [3] . 
In these two concentrated samples, the hysteresis cycles 
also present more complex features and are less repro- 
ducible so that in the following we shall mostly focus on 
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FIG. 2: (a) Amplitude £o of the surface waves versus driving 
acceleration a in suspensions of fd viruses vibrated at / = 
100 Hz for three different concentrations: (a) c = 5.8, (b) 
11.3, and (c) 13.6 mg/mL. The black line in (a) is the best fit 
by Eq. © with a = 42.9 m.s" 2 , n s = 10~ 3 Pa.s, r/i = 5.9r) s , 
?3 = 1.17 mm, and is = 1.3 mm. Dotted lines indicate a" p 
(see text). 
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FIG. 3: (a) Modulus of the complex viscosity \rj*\ vs. ampli- 
tude 70 of the oscillatory shear rate for the isotropic suspen- 
sion at c = 5.8 mg/mL. Colors refer to different oscillation 
pulsations: uj = 1 (•), u> = 3 (•), u> = 10 (•), uj = 30 (•), 
and uj — 100 rad.s -1 (•). The corresponding continuous lines 
are the best fits by Eq. © with 7i s = 10~ 3 Pa.s. (b) Char- 
acteristic time ti and (c) viscosity ni extracted from the fits 
by Eq. @ versus uj (o) and their extrapolations at 50 Hz (□) 
from quadratic fits in logarithmic scales (red lines). 



modelling the behaviour of the dilute sample. In any 
case, these effects are characteristic of a subcritical in- 
stability and can be attributed to the shear-thinning be- 
haviour of our rodlike colloids. A simple argument relat- 
ing shear-thinning to subcriticality is that a small pertur- 
bation of the surface may induce a shear rate that is large 
enough to significantly decrease the local viscosity, giving 
rise to a positive feedback on the perturbation and hence 
leading to finite wave amplitude at instability onset. As 
recalled above in the introduction, such an argument was 
the basis for a recent model by Epstein and Deegan [l2[ 
who used piecewise linear fits of the viscosity vs. shear 
rate curve to account for high-amplitude strip waves in 
dilute wormlike micelle solutions. It was also shown that 
shear-thinning leads to a hysteretic, discontinuous tran- 
sition. 



However, the model of Ref. [12[ was based on steady- 
shear rheology whereas the flow induced by Faraday 
waves is oscillatory. Therefore, in the present study, 
we propose a simple alternative approach based on rhe- 
ological measurements under large amplitude oscillatory 
shear. Figure [3] shows data recorded on the dilute sample 
using a standard rheometer (TA Instruments AR2000N). 
It is seen in Fig. Eta) that, whatever the oscillation pul- 
sation uj, the modulus of the complex viscosity if as a 
function of the amplitude of the oscillatory shear rate 70 
is well fitted by the Cross model 



21|: 



h*l(7o) 



1 + (ioriY 



(2) 



where t\ is a characteristic time, m an exponent, and 
•q s the solvent viscosity taken equal to the water viscos- 
ity (t] s = 10~ 3 Pa.s). The coefficient rji can be inter- 
preted as 771 = r/o — ?7sj where r/o is the modulus of the 
zero-shear complex viscosity at uj. For all fits, we get 
to = 1.5 ± 0.1 and the dependence of t\ and r)i on uj 
is shown in Fig. [3jb) and (c) respectively. Since most 
classical rheometers are limited to oscillation frequencies 
below about 10 Hz, we use quadratic fits in logarithmic 
scales to extrapolate both t\ and 771 to 50 Hz, i.e. to 
half the driving frequency used in the Faraday experi- 
ment [see red lines and green symbols in Fig. [3jb) and 
(c)]. Finally, in order to keep the analytical expressions 
tractable while including the above shear-thinning rhe- 
ology, we shall thereafter replace the exponent to = 1.5 
by to = 2 in Eq. $Z§. We do not expect this approxi- 
mation to change the results qualitatively but it should 
be noted that it sharpens the decrease of the viscosity 
around the characteristic shear rate 71 = I/77 = 235 s _1 
extrapolated at 50 Hz. 

Coming back to Faraday waves, the oscillating shear 



rate at the surface reads 7 = 



jLUt/2 



/2, where k c 



is the critical wave number of the surface pattern (ob- 
tained by considering half a mode in the y direction 22| : 
k c = \Jk' 2 + n 2 /I 2 . For the dilute sample, amplitude 
measurements give £0 ~ 6.10 -4 m and k c ~ 1380 m -1 at 
/ = 100 Hz, yielding | -y | ~ 260 s -1 . This shear rate 
is comparable to the characteristic shear rate 71 [see 
Fig. [3ja)] so that a strong effect of Faraday waves on 
the fluid viscosity is expected. At this stage, a model 
based on the numerical integration of the Mathieu equa- 
tion could be developed following Ref. [l2j ■ Here, how- 
ever, we propose to stick with analytical expressions by 
considering our shear-thinning fluid as an effective New- 
tonian fluid. More precisely, provided the fluid layer is 
deep enough (k c h 3> 1) and the viscosity is low enough 
(2rjk 2 /uj <S 1), which is indeed the case in our exper- 
iments where k c h ~ 13 and 2r)k 2 /u> ~ 0.06, the local 
shear rate j(z,t) — jo(z)eJ LJt ^ 2 is known to decrease ex- 
ponentially with the depth z: 



iok c u> 



(3) 
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where —h < z < by convention Q. Using Eq. ([2]) with 
m = 2 yields the z-dependence of the viscosity: 



\V*\(Z) 
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~ 7 T 



(4) 



Next, we define r] e s, the viscosity of a Newtonian fluid 
which dissipates the same energy as our colloidal suspen- 
sion for a given surface wave amplitude £o> as: 



f/eff 



Loo v{z)i{z) 2 dz 

J°oo7(^) 2 dz 



(5) 



where i c /i > 1 was used to extend the integration from 
—h to — oo. Introducing the characteristic amplitude 
£ c = 27i/(wfc c ), Eqs. ((H) and © readily lead to: 
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(6) 



Finally, for Newtonian fluids, it is well known that the 
amplitude of the surface wave should evolve with time 
as [231 : 
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where a c is the critical acceleration, (a — a c )/(a c r) is 
the growth rate, and ^3 and 1$ are characteristic lengths 
linked to the dissipation process. For Newtonian fluids 
of low viscosity, it was also shown that the critical ac- 
celeration is proportional to the viscosity [22I ]: a c oc 7y c ff. 
Therefore, at equilibrium (d/dt = 0), combining Eqs. ([5]) 
and j7]) yields the relationship between a and £o : 
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where a is the critical acceleration of the fluid at rest 
(i.e. for £0 = 0) and corresponds to the acceleration 
a" p measured experimentally. When £ — > 0, the second- 
order expansion of the previous expression in £g reads: 
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so that the instability is supercritical for Z3 



(9) 
< 



£ c a/2(1 + rjs/rji) and subcritical otherwise. In other 
words, when £ c is small enough or equivalently when 
shear-thinning sets in at a low enough shear rate 71 , sub- 
criticality and hysteresis are expected in the Faraday in- 
stability. Figure [2ja) shows that Eq. ([8]) provides a good 
fit of the experimental a(£o) data with only two free pa- 
rameters, ^3 and I5, since rj s is known, ao is measured 
independently as a" p , and 771 is extracted from Fig. EJc). 
Still, the predicted curve goes to zero for an acceleration 
smaller than the experimental aJ? own . We believe that 
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FIG. 4: From top to bottom: minimum of the transmitted 
intensity min(J), time-averaged intensity (J), and amplitude 
Ii of the first harmonic of the intensity recorded in suspen- 
sions of fd viruses vibrated at / = 100 Hz for three different 
concentrations: c = 5.8, 11.3, and 13.6 mg/mL from left to 
right. The white dotted lines indicate the bottom of the cell. 
All movies were taken at a — a" p . Intensities are normalized 
by the maximum intensity of a given movie. ^ 



the main reasons for this discrepancy are the very simple 
form of Eqs. d3| and ([5]) along with possible dissipation 
at the walls induced by the presence of a meniscus, which 
is not taken into account in our simple model. We shall 
come back to the two more concentrated samples in the 
next section. 



BIREFRINGENCE MEASUREMENTS 

The above model provides a link between hysteresis 
and shear-thinning induced by the alignment of the rod- 
like colloids. This alignment is directly confirmed ex- 
perimentally in Fig. [4] that reports the analysis of the 
transmitted intensity field I(x, z, t) recorded over twenty 
driving periods with crossed polarizers:' the minimum 
intensity and the time-averaged intensity are presented 
in the first two rows while the last row shows the am- 
plitude 1\ of the first harmonic of I(x,z,t), i.e. of the 
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component of / that oscillates at w. In all cases, bright 
spots are observed near the surface that are evidence of 
a localized alignment of the fd viruses. More precisely, as 
long as the system remains far from a phase transition, 
the transmitted intensit y fi eld can be interpreted using 
the "stress optical rule" [23| according to which the bire- 
fringence An is proportional to the shear stress tr, so that 
/ oc An 2 oc a 2 . Therefore, at least for the more dilute 
sample which is far from the I-N transition, the birefrin- 
gence intensity should oscillate at u) and the amplitude 
of the first harmonic should read: 



I\{x, z) — KS\n 2 {kx) |n*(2)| 2 \{z) . 



(10) 



where k is a constant that depends on the exact geom- 
etry of the experiment, the optical set up, and the op- 
tical properties of the sample. Figure [5ja) shows that 
the transmitted intensity indeed oscillates at u. More- 
over, for the dilute sample, I\{x, z) can be predicted 
from the model developed above by inserting Eqs. ([3]) 
and (J4j> into Eq. (fT0|) . We restrain ourselves to the 
z-dependence of the transmitted intensity by calculat- 
ing the average over x of the temporal standard devi- 
ation of the transmitted intensity I(x,z,t): SI(z) =< 
y/< (I(x,z,t)- < I(x,z,t) > t ) 2 >t >x (where <>. i de- 
notes the average over i), which should be proportional 
to |?7*(z)| 2 7q(z). As seen in Fig.EJb), this approach pro- 
vides a good description of 51{z) with the constant k 
as the only free parameter. Also, when plotted versus 
z — ln(£o/£c)/fcc ; all curves 51 (z) measured for different 
vibration amplitudes collapse on the same curve. Inci- 
dentally, this also proves that k does not depend on £ - 
More importantly, the results shown in Figs. HJleft) and 
[5] for the more dilute sample are consistent with Eq. (fTU)) 
and accredit a scenario of shear-induced alignment in 
which the rodlike colloids simply orient periodically ac- 
cording to the shear-rate field generated by the surface 
wave pattern. 

For samples close to or above the I-N transition, how- 
ever, the stress optical rule breaks down as the or- 
der parameter saturates and the above approach is no 
longer justified. Still, interesting new effects show up 
in Fig. @|middle) for the sample at c = 11.3 mg/mL. 
There, it is seen that the minimum intensity is nonzero 
over a thin layer below the surface covering almost the 
whole width of the sample. This layer also corresponds 
to maxima in (I) while it does not show in I\. This 
means that the surface waves have generated a perma- 
nent birefringence pattern below the surface with the 
same wavelength. t In other words, just below the I-N 
transition, the sample keeps the imprint of the surface 
wave pattern. The fact that the continuous component 
of the transmitted intensity becomes larger for larger con- 
centrations is also reflected in the smaller oscillations of 
the relative amplitude reported in Fig. [5l[a): in the more 
concentrated samples, the oscillating component repre- 
sents less than 5 % of the average intensity, about three 
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FIG. 5: (a) Normalized transmitted intensity I /(I) — 1 versus 
time t recorded at an intensity antinode in suspensions of fd 
viruses vibrated at / = 100 Hz for three different concentra- 
tions: c = 5.8 (red), 11.3 (green), and 13.6 mg/mL (blue), and 
a = 1.1 a" p . (b) 5I(z) versus z — ln(£o / £c) / 'k c recorded in the 
more dilute suspension (c = 5.8 mg/mL) at 100 Hz and for 
different vibration amplitudes (see text). The grey line rep- 
resents Eq. (| 10p where jo(z) is given by Eq. @ and |?7*|(.z) 
by Eq. Q. Except for the prefactor k, all the parameters are 
inferred from the rheological measurements of Fig. [3] 



times smaller than for the dilute sample. 

As for the most concentrated sample at c = 
13.6 mg/mL, it is easily checked from (I) in 
Fig. BJmiddle) that this nematic phase presents some 
significant birefringence everywhere in the sample. More- 
over, although the effect is less spectacular than for the 
sample at c = 11.3 mg/mL due to this natural birefrin- 
gence, a continuous pattern is also generated that reflects 
the surface wave pattern. We interpret the build-up of a 
continuous birefringence field as a trade-off between rota- 
tional relaxation of the colloids and shear-induced align- 
ment. In steady state, the rodlike viruses align locally 
due to shear, but close to the I N transition, their large 
relaxation time prevents them from reorienting between 
two oscillations so that they remain aligned. This quali- 
tative argument, which is developed more formally below, 
may also explain why the hysteresis cycles of the concen- 
trated suspensions could not be properly accounted for 
[see Fig. [2|b) and (c)]: since it takes longer for the fluid 
to relax after alignment, our experimental determination 
of the critical acceleration may not be accurate enough 
leading to some irreproducibility in a c . 



NUMERICAL ANALYSIS 

A basic understanding of the virus alignment induced 
by Faraday waves can be obtained through numerical 
analysis following the calculation by Dhont and Briels 
for rig id cylindrical particles of length L and diameter 
D I n short, this approach yields the temporal 

evolution of the orientational order parameter tensor S 
based on the shear-rate field. One result is that, in the 
isotropic phase (L<fi/D < 4), the linear complex viscosity 



6 



r]*(u)) follows the Oldroyd-B model, where the charac- 
teristic variables depend on the solvent viscosity r) s , the 
volume fraction of rodlikc particles </>, their aspect ratio 
L/D, and their rotational diffusion coefficient at infinite 
dilution D r , according to: 



L(])/D=2 Li|>/D=4 



rf{u) = rj s + 



1 + iuiT 

1 

3rj s a 



T = \6D r 1 - 



L(f> 
~5D 



with a = — | — — 



5 \3D J ln(L/D) 



(11) 
(12) 
(13) 
(14) 

(15) 



It should be noted that the above expression is only 
valid in the isotropic state and that extra care should 
be taken in the nematic state by fitting the calculated 
viscosity. Under this model, the sample is isotropic for 
L<f>/D < 40/9 and nematic for L4>/D > 40/9. More- 
over no direct comparison can be made between rheo- 
logical results from this model and experimental ones as 
this model neglects the flexibility of the rodlike colloids 
which becomes important at high frequency due to Rouse 
modes, while at low frequency the particle diffusion is 
neglected 25]. Still this calculation has proven useful 
to predict the alignment of rigid rods under shear (20j . 
which we investigate here. 

In the limit of small amplitude (linear approximation) 
in an horizontally unbounded cell, Ref. Q gives the nu- 
merical expression of the two-dimensional velocity field 
v(a:, z) regime knowing the complex viscosity of the sam- 
ple. Using Eq. (|TT1) . we first calculate the v(x,z) in the 
linear regime. We then deduce the shear rate field from 
which we calculate the local order parameter tensor S 
following Refs. 19|, 24 1 with only one additional input, 
namely the amplitude £o of the surface wave. 

We also compute the birefringence intensity I from the 
order parameter S: 



I = h 



sin(2(0 - 6>)) 2 (1 - cos(27r/Ai//A)) 



(16) 



with Av = — [ I ' Aa ^(s xx -s zz y+Asun) 

18 n 



, n 1 ( 2 Sxz 

and — — arctan — 

2 \ S X x 



G18) 



where the input parameters are the orientation of the 
pair of crossed polarizers /3, the sample thickness I, the 
incident light wavelength A = 575 nm chosen as yellow 
(although we used white light experimentally), the nor- 
malized concentration v = 40/ '(irD 2 L), the mean refrac- 
tive index approximated as the water index n = 1.33, 
and the polarizability anisotropy Aa = 1.5.10 -25 m -3 
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FIG. 6: From top to bottom: minimum of the order parameter 
minus its isotropic value A — 1/3, time-averaged birefringence 
intensity (/), and amplitude Ii of the first harmonic of the 
intensity calculated for three concentrations approaching the 
I-N transition L<f>/D — 2, 4, and 4.44, and a concentration in 
the nematic phase L<j>/D — 4.67. 



deduced from Ref. [26| . For the calculations of Fig. [51 we 
take /3 = 45°, I = 6 mm, and £o = 600 /.tm consistently 
with the experimental values. 

From these numerical calculations, we extract the min- 
imum of the order parameter A, the mean intensity (/), 
and the amplitude of the first harmonic I\ of the intensity 
field (respectively shown from top to bottom in Fig. 
for three isotropic suspensions (L/Dcf> = 2, 4, and 4.44) 
and one in the nematic phase {L/Dcj) =4.67). The mini- 
mum of the order parameter \{x, z) (first row of Fig. [5]) 
is calculated as the minimum during an oscillation of the 
maximum eigenvalue of S(x, z,t) and should be equal to 
1/3 for isotropic systems and greater than 1/3 in the case 
of an out-of-equilibrium I-N transition. Close to the sur- 
face, a band is present where A > 1/3 with a periodicity 
2k. This is a clear indication of a phase transition in- 
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duced by the surface waves and is similar to the bright 
zones of the minimum of intensity seen in Fig. |4j The 
various features of the experimental birefringence pattern 
are qualitatively well captured. Especially, the width of 
the aligned band increases when approaching the I-N 
transition. This is a signature of the slowing down of 
the system at the I-N transition: as the relaxation time 
increases, the critical shear rate needed to align the sus- 
pension decreases. 

Although the qualitative features of Fig. @] are well re- 
produced by the numerical calculations, significant dis- 
crepancies are observed, most importantly the presence 
of higher order of birefringence (2irlAv/A > it) result- 
ing in a secondary maximum in the computed intensity 
field that is absent in the experiment. This intensity 
peak marks the end of higher order birefringence. At the 
position where this maximum occurs, one has A = 1/3. 
Therefore, at this point, the sample is not permanently 
aligned. Since all the pref actors entering Ais are well 
known [see Eq. (j!7[) ] , we conclude that the order parame- 
ter S is different in the experiment and in the calculation. 
Many reasons can lead to this discrepancy. As already 
mentioned, the rheology deduced from Refs. 3, 24 1 dif- 
fers from the experimental rheological behaviour which 
leads to different velocity profiles. The calculation also 
supposes that the deformations remain small (linear hy- 
pothesis). As such, we should have £ok •C 1 which is 
not the case in our experiments where ^ok c ~ 1. Thus, 
in order to reach quantitative agreement, nonlinear terms 
related to the finite amplitude of the surface waves should 
be taken into account while calculating the velocity field. 
Finally, the presence of a meniscus at the edge of the 
cell can damp the velocity. Still, the basic calculations 
presented in Fig. [6] strongly support our interpretation 
of the experimental observations in terms of a localized 
isotropic-to-nematic transition triggered by the Faraday 
instability. 



CONCLUSION 

The present study has unveiled a strong effect of the 
Faraday instability on suspensions of rodlike colloidal 
particles. By measuring the surface height profile, we 
have shown that the transition to parametrically excited 
surface waves displays discontinuous, hysteretic features. 
We first linked this subcritical behaviour to the shear- 
thinning properties of our colloidal suspensions thanks 
to a phenomenological model based on rheological data 
under large amplitude oscillatory shear. We have fur- 
ther provided evidence that Faraday waves induce local 
nematic ordering of the rodlike colloids through bire- 
fringence measurements. While local alignment simply 
follows the surface oscillations for dilute, isotropic sus- 
pensions, permanent nematic patches are generated by 
surface wave in samples close to the isotropic-to-nematic 



transition. Above the I-N transition nematic domains 
align in the flow direction. Such a strong coupling be- 
tween the fluid microstructure and a hydrodynamic insta- 
bility was confirmed by numerical computations. These 
results clearly differ from what was already observed in 
semidilute wormlike micelles Q where the patterning of 
the birefringence intensity came from stretching of the 
micellar network without any phase transition. 

Although more experiments are needed to fully un- 
derstand the behaviour of rodlike colloids in complex 
shear fields, our results show that localized isotropic-to- 
nematic phase transitions can be induced in temporally 
stable structures with the same spatial organization as 
the surface instability pattern. This may open new paths 
of research, both experimental and theoretical, as well 
as potential applications, in which a complex system is 
structured at a mechanically or hydrodynamically con- 
trolled length scale intermediate between the mesoscopic 
size of the microstructure and the macroscopic size of the 
container. 

The experiments were performed while the authors 
were at Centre de Recherche Paul Pascal (CRPP). We 
wish to thank the "Cellule Instrumentation" of CRPP 
for technical advice and design of the experiment as well 
as E. Grelet, S. Lerouge, and J. Dhont for enlightening 
discussions. 
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